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We compute the masses of the singly and doubly charmed baryons in full QCD using the 
relativistic Fermilab action for the charm quark. For the light quarks we use domain-wall 
fermions in the valence sector and improved Kogut-Susskind sea quarks. We use the low- 
lying charmonium spectrum to tune our heavy-quark action and as a guide to understanding 
the discretization errors associated with the heavy quark. Our results are in good agreement 
with experiment within our systematics. For the S^c, we find the isospin-averaged mass to 
be Mh„ = 3665 ± 17 ± 14 MeV; the three given uncertainties are statistical, systematic 
and an estimate of lattice discretization errors, respectively. In addition, we predict the 
mass splitting of the (isospin-averaged) spin-1/2 flee with the Sec to be Afji^c ~ ^'k^c — 
98 ± 9 ± 22 ± 13 MeV (in this mass splitting, the leading discretization errors are also 
suppressed by SU{3) symmetry). Combining this splitting with our determination of Mh^^ 
leads to our prediction of the spin-1/2 ilcc mass, Ma„ = 3763 ± 19 ± 26 t^l MeV. 



I. INTRODUCTION 



Experimental and theoretical studies of charmed and bottom hadrons have been the focus of vig- 
orous research over the last several years [H EJ [Sj |3j . In particular, singly and doubly heavy baryon 
spectroscopy has received significant attention, mainly due to the recent experimental discoveries 
of both new charmed (SELEX) [5l [6] and bottom baryons by DO [4 and CDF [7 . In addition to 
these discoveries, there are still many states of heavy and doubly heavy baryons remaining to be 
discovered. The new Beijing Spectrometer (BES-III), a detector at the recently upgraded Beijing 
Electron Positron Collider (BEPCII), has great potential for accumulating large numbers of events 
to help us understand more about charmed hadrons. The antiProton ANnihilation at DArmstadt 
(PANDA) experiment, a GSI future project, and the LHCb are also expected to provide new re- 
sults to help experimentally map out the heavy-baryon sector. For these reasons, lattice quantum 
chromodynamic (QCD) calculations of the spectrum of heavy baryons are now very timely and will 
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play a significant role in providing theoretical first-principles input to the experimental program. 

Lattice QCD is now a mature field capable of providing accurate results that can be directly 
compared to experiment, with calculations in the light-quark sector being well established. Al- 
though the study of heavy quarks requires careful treatment of discretization errors, significant 
advances have been made in this sector as well. Lattice heavy quarks have ©((mga)") errors, 
where mg is the mass of the heavy quark and a is the lattice spacing. Lattice spacings for typ- 
ical, currently accessible dynamical ensembles are still too coarse (a~^ « 2 GeV) to make such 
systematic errors small. To assert better control over the discretization errors for heavy quarks on 
the lattice, several heavy-quark approaches have proven useful. For example, non-relativistic QCD 
(NRQCD) [8:, which is an expansion of the lattice quark action in powers of is commonly 

applied to bottom quarks. However, the charm-quark mass is not heavy enough to justify the use of 
NRQCD. Relativistic heavy-quark actions [H [101 El US] systematically remove ©((mga)") terms 
and are better suited to charm-quark calculations. Recent updates on the state of heavy-quark 
physics on the lattice can be found in several reviews |13| [T^ [T5l [TBI El HH] ^-nd references therein. 

Up to now, there have been a few lattice charmed-baryon calculations using the quenched ap- 
proximation. In some cases an C'(a)-improved light-quark action is used on isotropic or anisotropic 
lattices with a single lattice spacing: Bowler et al. p2j used a tree-level clover action for both 
light and heavy quarks to calculate the singly charmed baryons spectrum of spin 1/2 and 3/2. 
Later, Flynn et al. [20] updated this project with nonperturbative clover action and extended 
the calculation to doubly charmed baryons. Chiu et al. [21; used a chiral fermion action for the 
charm quarks and calculated both the positive and negative parity spectrum for singly and doubly 
charmed baryons. Such calculations using light-quark actions to simulate heavy quarks introduce 
large systematic errors proportional to (amg)^, which must be carefully addressed. One calcula- 
tion has used a higher-order improved fermion action: Lewis et al. |22J performed a calculation 
on both doubly and singly charmed baryons using D234-type fermion action (which would leave a 
leading error of 0{a'^)) for both light and heavy quarks but on a coarse anisotropic ensemble (with 
anisotropy = 2). Finally, heavy-quark effective theory was applied to charm calculation: Mathur 
et al. [22 continued to use anisotropic lattices, adding two more lattice spacings, but changed the 
heavy-quark action to NRQCD, which reduces the lattice-spacing discretization effects. For all 
of these calculations, the quenched approximation remains a significant source of systematic error 
that is difficult to estimate. 

Given the progress on the experimental side, it is time to revisit these charmed baryon calcu- 
lations using dynamical gauge ensembles and improve the calculations with the current available 
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computational resources. Although more dynamical ensembles are available these days, not many 
charmed baryon calculations have been published so far, only a few proceedings |24l l25l |26] . 

In this work, we extend our previous calculation [SB] to higher statistics and compute the ground- 
state spectrum of the spin-1/2 singly and doubly charmed baryons. We use the Fermilab action [9] 
for the charm quarks and domain-wall fermions for the light valence quarks on gauge configurations 
with 2-|-l-flavor Kogut-Susskind fermions and a range of quark masses resulting in pion masses as 
light as 290 MeV. We nonperturbatively tune the fermion anisotropy and two input bare masses 
for charm quarks, setting the remaining parameters to tree-level tadpole improved coefficients. 
Our results are extrapolated to the physical light-quark masses using both heavy-hadron chiral 
perturbation theory (HHxPT) as well as HHxPT-inspired polynomial extrapolations. 

II. LATTICE FORMULATION 
A. Light-Quark Action 

In this work we employ the "coarse" (a ~ 0.125 fm) gauge configurations generated by the MILC 
Collaboration [27] using the one- loop tadpole-improved gauge action p8], where both 0{a'^) and 
O^g^a?) errors are removed. For the fermions in the vacuum, the asqtad-improved Kogut-Susskind 
action [29l EO] EU ISl [Ml Ej is used. This is the Naik action [35^ (©(a^) improved Kogut-Susskind 
action) with smeared links for the one-link terms so that couplings to gluons with any of their 
momentum components equal to vr/a are set to zero. 

For the valence light quarks (up, down and strange) we use the five-dimensional Shamir |36| I37| 
domain- wall fermion propagators [38j calculated by the NPLQCD Collaboration [39] . The domain- 
wall fermion action introduces a fifth dimension of extent L5 and a mass parameter M5; in our 
case the values L5 = 16 and = 1.7 were chosen. The physical quark fields, q{x,t), reside on 
the 4-dimensional boundaries of the fifth coordinate. The left and right chiral components are 
separated on the corresponding boundaries, resulting in an action with chiral symmetry at finite 
lattice spacing as L5 ^ cxd. We use hypercubic-smeared gauge links |40| E] H2| H3] to minimize 
the residual chiral symmetry breaking, and the bare quark-mass parameter (am)^'"^ is introduced 
as a direct coupling of the boundary chiral components. 

The calculation we have performed, because the valence and sea quark actions are different, 
is inherently partially quenched and therefore violates unitarity. Unlike conventional partially 
quenched calculations, to restore unitarity, one must take the continuum limit in addition to 



4 



Ensemble 


P 


ami 


arris 






-^cfgs 


props 


mOOr 


6.76 


0.007 


0.050 


0.0081 


0.081 


461 


2766 


mOlO 


6.76 


0.010 


0.050 


0.0138 


0.081 


636 


3816 


m020 


6.79 


0.020 


0.050 


0.0313 


0.081 


480 


1920 


m030 


6.81 


0.030 


0.050 


0.0478 


0.081 


563 


1689 



TABLE I: The parameters of the configurations and domain-wall propagators used in this work. The 
subscript I denotes light quark, and s denotes the strange quark. The superscript "dwf" denotes domain- 
wall fermion. 

tuning the valence and sea quark masses to be degenerate. This process is aided with the use of 
mixed-action chiral perturbation theory |44| H5l H6| IT7| SHI HH]. Given the situation, there is an 
ambiguity in the choice of the valence Hght-quark masses. One appeahng choice is to tune the 
masses such that the valence pion mass is degenerate with one of the staggered pion masses. In the 
continuum hmit, the Nf = 2 staggered action has an SU{8)l 5'C/(8)/j (8) U{l)v chiral symmetry 
due to the four- fold taste degeneracy of each flavor, and each pion has 15 degenerate partners. At 
finite lattice spacing this symmetry is broken and the taste multiplets are no longer degenerate, but 
have sphttings that are 0{ala^) [291 EQl EH Si ED] . The propagators used in this work were tuned 
to give valence pions that match the Goldstone Kogut-Susskind pion. This is the only pion that 
becomes massless in the chiral limit at finite lattice spacing. As a result of this choice, the valence 
pions are as light as possible, while being tuned to one of the staggered pion masses, providing 
better convergence in the needed to extrapolate the lattice results to the physical quark-mass 
point. This set of parameters, listed in Table |l| was first used by LHPC |5H [52] and recently to 
compute the spectroscopy hadrons composed of up, down and strange quarks [53]. 

B. Heavy- Quark Action 

For the charm quark we use the Fermilab action [9], which controls discretization errors of 
©((amg)"). Following the Symanzik improvement |53j, an effective continuum action is con- 
structed using operators that are invariant under discrete rotations, parity-reversal and charge- 
conjugation transformations, representing the long-distance limit of our lattice theory, including 
leading finite-a errors. Using only the Dirac operator and the gluon field tensor (and distinguishing 
between the time and space components of each), we enumerate seven operators with dimension 
up to five. By applying the isospectral transformations [55], the redundant operators are identified 
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and their coefficients are set to appropriate convenient values. The lattice action then takes the 
form 



S — Sq + Sb -\- Se , 



(1) 



with 



5o = ^Q(a 



mo + (70 Vo - -Ao) + (7^V^ - - A 



Q{x), 



Sb = -^CB^Qix) ^CTijFij Q{x) , 



i<j 



Se = -^CE^Q{x) ^^<ToiFoi^ Q{x) , 



(2) 
(3) 
(4) 



where a is the lattice spacing, Vo and Vj are first-order lattice derivatives in the time and space 
directions, Aq and Aj are second-order lattice derivatives, and F^^, is the gauge field strength 
tensor. The spectrum of heavy-quark bound states can be determined accurately through \p\a and 
{amo)'"' for arbitrary exponent n by using a lattice action containing mo, cb and c^;, which are 
functions of amq. 

The coefficients cb and ce are different due to the broken space-time interchange symme- 
try, which can be computed in perturbation theory by requiring elimination of the heavy-quark 
discretization errors at a given order in the strong coupling constant as- We use the tree-level 
tadpole-improved results obtained by using field transformation (as in Ref. |55j): 



cb 



CE 



(5) 







where uq is the tadpole factor 



no 



1/4 



(6) 



and Up is the product of gauge links around the fundamental lattice plaquette p. The remaining 
two parameters m-o and u are determined nonperturbatively. The bare charm-quark mass uiq is 
tuned so that the experimentally observed spin average of the J /'^ and r]c masses 



avg 



(7) 



is reproduced; see Sec. IV B for further details. The value of v must be tuned to restore the 
dispersion relation Ef^ = m| -|- (P'p^ such that = 1. Since the values of v and mo are coupled, 
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Ensemble 


c2 


Vc 




D 


Ds 


m007 


0.991(4) 


0.985(5) 


1.021(15) 


1.018(9) 


mOlO 


0.989(3) 


0.958(3) 


1.016(10) 


0.992(6) 


m020 


0.997(4) 


0.993(5) 


1.019(20) 


1.004(14) 


m030 


0.963(5) 


0.947(6) 


1.029(12) 


1.015(10) 



TABLE II: Speed of light for charmed mesons. 

one needs to iterate the tuning process in order to achieve a consistent pair of values. To do this, 
we calculate the single-particle energy of rjc, J/^, Dg and D at the six lowest momenta (with unit 
ofa-i): ^(0,0,0), ^(1,0,0), ^(1, 1, 0), ^(1, 1, 1), ^(2, 0, 0), ^(2, 1, 0). For each ensemble, the 
energy levels are calculated at two charm-quark masses (denoted mi = 0.2034 and m2 = 0.2100) 
and extrapolated to the physical charm-quark mass (as described below). The values of are 
obtained by fitting the extrapolated energy levels to the dispersion relation. We tune u using the 
dispersion relation of ijc- As one can see from Table the dispersion relations for either the 
charmonium J/^* or the charm-light mesons (D and Dg) are generally consistent with = 1 to 
within 1-2%. 



III. CHARMED HADRON SPECTRUM: NUMERICAL RESULTS 



The interpolating operators we use for the J = 1/2 singly and doubly charmed baryons are 



A. 



JjkfJT 



JjkfJT 



1 



ijk 



quCj,Qi)q'g+{qfCj5Qi)q'u 



V2 

e'^HqfC^5Qi)q's, 
e'^\CrIC^^qi)Ql, 
e''HQfCj,qi)Q 



k 



(8) 



where q^^d ^-re the up and down quark fields, qs is strange quark field and Qc is charm quark field. 
Using these interpolating fields, we construct the two-point functions 
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where Oh is an interpolating operator of the hadron h. The correlation functions are calculated 
with gauge-invariant Gaussian-smeared sources and point sinks. The smearing parameters were 
optimized so that excited-state contamination to the correlators is minimized. The domain-wall 
valence propagators were computed with Dirichlet boundary conditions in the time direction, re- 
ducing the original lattices to half their temporal size. Similar to baryons, the signal for the 
charmed correlation functions quickly drops, and thus we do not expect the temporal reduction to 
reduce the number of useful time points for our analysis. The sources were located away from the 
Dirichlet boundary to minimize contamination from the boundary effects. In order to enhance our 
statistical precision, several valence propagators are taken from each configuration with varying 
source location. The resulting correlation functions are then source averaged on each configuration 
to produce one correlator per configuration for each interpolating operator. The masses of the 
hadrons are obtained by fitting the correlation functions to a single exponential 

Chit) = ^e-^o* (9) 

in a region where the effective mass is observed to exhibit a plateau. The fitting range is varied to 



estimate the systematics from the choice of fitting window, as indicated in Tables |III| and IV 



In Table III we summarize the resulting baryon masses as well as the corresponding time ranges. 



The first uncertainty is statistical and the second is a fitting systematic. For most fits, the resulting 
per degree of freedom is about one. In Figure [l] we display representative effective mass plots 
and their fitted masses for both good and poor fits. The results from charmonium are shown in 
Table HVl 



IV. HEAVY- AND LIGHT-QUARK MASS EXTRAPOLATION 

In order to make contact with experiment, we must extrapolate our results to infinite volume, 
continuum limit and to the physical value of the light- and heavy-quark masses. Optimally, the 
extrapolations can be performed in terms of dimensionless ratios of observable quantities, so as 
to minimize contamination from a particular scale-setting method. In this work, we have chosen 
to scale our masses by the calculated value of the pion decay constant on each ensemble, forming 
the dimensionless ratios Mh/fn, where Mh is the mass of a given hadron. We take the values of 

(and rTT-Tr) from Ref. |53j; they are collected in Table [V} As can be seen, afj^ varies by ~ 15% 
over the range of pion masses used in this work, adding additional chiral curvature. However, 
the light-quark mass dependence of is well understood |56| l57] . and so this variation can be 
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Hadron 


mo m007 mOlO m020 m030 




mi 2.3578(18)(8)[8-18] 2.3620(14)(9)[10-18] 2.3456(33)(17)[12-18] 2.3333(23)(6)[1H8] 
7712 2.3663(18)(8)[8-18] 2.3705(14)(9)[10-18] 2.3542(33)(16)[12-18] 2.3419(23)(7)[11-18] 


^cc 


7771 2.3018(27)(0)[7-13] 2.3120(23)(23)[9-17] 2.3087(33)(3)[8-18] 2.3056(28)(33)[11-18] 

7772 2.3104(27) (0) [7-13] 2.3205(23) (23) [9-17] 2.3173(33)(3)[8-18] 2.3142(28)(33)[11-18] 




mi 1.7216(24)(1)[9-15] 1.7240(24)(5)[12-18] 1.7101(52)(77)[12-16] 1.7160(39)(13)[12-18] 
m2 1.7261(24)(1)[9-15] 1.7285(24)(5)[12-18] 1.7146(52)(76)[12-16] 1.7205(39)(13)[12-18] 


'— 'c 


7771 1.6754(26)(32)[6-18] 1.6799(29)(43)[9-16] 1.6875(52) (57) [9-16] 1.6881(43)(2)[11-18] 

7772 1.6799(26)(32)[6-18] 1.6844(29)(43)[9-16] 1.6920(52) (58) [9-16] 1.6927(43)(2)[11-18] 




7771 1.6076(82)(86) [12-18] 1.6078(48)(54)[12-18] 1.6167(40) (9) [8-18] 1.6120(41)(47)[12-17] 

7772 1.6121(82)(87) [12-18] 1.6123(48)(55)[12-18] 1.6211(40)(9)[8-18] 1.6163(41)(48)[12-17] 




7771 1.6157(50)(38)[7-17] 1.6252(55(0)) [9-15] 1. 6446(56) (0) [8-16] 1.6661(43)(70)[10-18] 

7772 1.6203(50)(38)[7-17] 1.6298(55)(0)[9-15] 1. 6491(56) (0) [8-16] 1.6706(43)(69)[10-18] 


Ac 


mi 1.4974(71) (47) [6-13] 1.523(16)(3)[12-18] 1.5571(55)(22)[8-18] 1.572(5)(18) [12-17] 
777,2 1.5018(71)(48)[6-13] 1.527(16)(3)[12-18] 1.5615(55)(22)[8-18] 1.577(5)(18) [12-17] 



TABLE III: Charmed baryon masses in lattice units with 2 values of mo (indicated as 7711 = 0.2034 and 
m2 = 0.2100) in Eq. ([2|. The first uncertainty is statistical and the second is systematic from the different 
choice of fitting ranges (presented in square brackets). 

accounted for. 

Ultimately, one would like to use lieavy-hadron chiral perturbation theory (HHxPT) \58\ \59\ \6U[ 
\6T\ [62l [63l to perform both the charm-quark mass extrapolation and the chiral extrapolation of 
the charmed hadron masses, allowing a lattice determination of not just the spectrum but also the 
low-energy constants entering the effective field theory. There are several reasons we cannot perform 
a thorough extrapolation in this manner. First, we only have results at four independent values 
of the light-quark mass, and at only one value of the strange mass. Second, in this work, we only 
have results for the J = 1/2 baryons, and a proper chiral extrapolation requires also the spectrum 
of J = 3/2 charmed baryons; the states are related by the heavy-quark symmetry, and therefore 
the mass splittings are small (similarly, the extrapolation of the heavy meson masses requires the 
J = 1 states as well as J = 0). Third, our calculation is mixed-action, thus requiring either a 
continuum extrapolation or the use of mixed-action xPT |44| H5l H6l ITTl SHj |39]. The mixed-action 
effective field theory can be trivially constructed from the partially quenched theories for heavy 
hadrons 165\ \66\ [67] by following the prescription in Ref. |48] . However, this work only utilizes one 
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Hadron 


m 


m007 


mOlO 


m020 


m030 


Vc 


mi 

7712 


1.8783(4)(0)[14-19] 
1.8866(4)(1)[14-19] 


1. 8804(3) (0) [12-19] 
1.8887(3)(1)[12-19] 


1.8687(4)(1)[12-19] 
1.8771(4)(1)[12-19] 


1.8598(3)(2)[8-15] 
1.8683(5)(0)[8-15] 


J/* 


mi 

1712 


1.9390(7)(0)[14-18] 
1.9470(7)(0)[14-18] 


1.9421 (4) (0) [10-19] 
1.9501(4)(1)[10-19] 


1.9296(6)(1)[12-19] 
1.9376(6)(1)[12-19] 


1.9198(6)(2)[11-19] 
1.9278(6)(3)[11-19] 


XcO 


mi 
m2 


2.1660(54)(21)[9-16] 
2.1741(54)(20)[9-16] 


2.1803(33)(6)[6-17] 
2.1883(35)(6)[6-17] 


2.1652(55)(50)[6-18] 
2.1733(55)(49)[6-18] 


2.1626(54)(2)[6-18] 
2.1705(54)(2)[6-18] 


Xcl 


mi 
m2 


2.2092(69) (24) [9-18] 
2.2171(69)(24)[9-18] 


2.2234(52)(35)[9-16] 
2.2312(52)(35)[9-16] 


2.2123(40)(8)[4-17] 
2.2199(40)(9)[4-17] 


2.2004(44) (25) [4-17] 
2.2081 (44) (25) [4-17] 


Xc2 


mi 

1712 


2.2224(64) (86) [6-18] 
2.2301(65)(85)[6-18] 


2.2386(32) (24) [4-18] 
2.2463(32) (25) [4-18] 


2.2205(45) (21) [4-17] 
2.2282(46)(19)[4-17] 


2.2151(63)(26)[5-18] 
2.2226(63)(25)[5-18] 



TABLE IV: Charmonium masses in lattice units with toi = 0.2034 and m2 = 0.2100. 



/3 

ensemble: 

ami 


6.76 
0.007 


6.76 
0.010 


6.79 
0.020 


6.81 
0.030 


am^^ 


0.1842 


0.2238 


0.3113 


0.3752 




0.0929 


0.0963 


0.1026 


0.1076 




1.983 


2.325 


3.035 


3.489 



TABLE V: Values of m^, and /^r calculated in Ref. [531. For all ensembles the staggered strange-quark mass 
is ams ~ 0.050 while the domain-wall strange-quark mass is am^"^ = 0.081. 

lattice spacing, and so one can not perform the full mixed-action analysis. With these caveats in 
mind, we proceed with our analysis. 



A. Scale setting with /^r 

The light-quark mass expansion of a heavy-hadron mass is given by^ 

(2) 

M. = Mo + ^?^ + ... (10) 

At this order, we are free to make the replacements /o f-w and 2Bmi m^, with corrections 
appearing at 0{m'^). The dots represent terms of higher order in the chiral expansion, with the 



^ Here we are presenting an SU{2) extrapolation formula with the operator normaUzation of Ref. |68) such that the 



(2) 

coefficient c\ is dimensionless. 
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:§ 1.60 ' 



1.85 



1.75 



10 

t/a 



10 

t/a 



2.38 - 

nj 2.34 - 

] 

2.32 - 
2.30 — 



Jj j rb 



§ 1.6 
CS 

1.5 



10 

t/a 



t/a 



FIG. 1: Sample efFective-mass plots and corresponding fits to the correlation functions. The smaller error 
bands are statistical and the larger error bands are statistical and systematic (determined by varying fit 
range) added in quadrature. 



first non-analytic (in the quark mass) corrections appearing as corrections which scale as 



As stated above, we are scaling our masses with f-,^ to form dimensionless ratios for extrapolation, 



Mh Mo , c 



(2) 



(11) 



When performing an extrapolation in this manner, it is important to realize we cannot approximate 
Mq/ as a constant, since the chiral corrections to are ©(m^) and thus are the same order 

.(2; 

normalization /q ^ 130 MeV) 



as the term with coefficient . Rather, the chiral expansion of is given by [56] (with the 



/tt — /o 



1 + 5f{mJU) 



+ ... 



+ 



(12) 



In this expression, we have made use of perturbation theory to replace all terms appearing at 
next-to-leading order with their (lattice) physical values. Similarly, we have rescaled the renormal- 
ization scale /i ^ /x/tt to express the chiral corrections as purely a function of m-,^/ f.,^. Again, the 
corrections to this rescaling first appear at next-to-next-to-leading order. In order to perform our 
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range 


m007-m010 


m007-m020 


m007-m030 


^4(M = /tt) 


0.0307(27) 


0.0293(6) 


0.0302(4) 



TABLE VI: Values of I4 needed for chiral extrapolations of Mh/ fn- The different values of I4 are determined 
through the different choices of fitting range, also listed. 




0.075 



m„|{2^^2nf„) 

(a) 



(b) 



FIG. 2: The (blue) filled circles represent the lattice data and the (red) star is the physical point, converted 
to lattice units using a^^ = 1588 MeV with a 2% error bar added for the scale setting. The error bands are 
the 68% confidence intervals in the resulting chiral extrapolation from the lightest two points (a) and a fit 
to all four lattice points (b). 



chiral extrapolations using Eq. (11), we must determine ^4, which captures the chiral corrections 



of /tt- The mixed-action formula for /^r is known [44J, but again, only useful if one has data for at 
least two lattice spacings. Since we currently only have results at one lattice spacing, we perform a 
continuum chiral extrapolation analysis of the c/tt in Table [V] The results are collected in Table [VTt 



The resulting extrapolations are plotted in Figure [2j In this figure, the (blue) filled circles are 
the lattice data, and the error bands represent the 68% confidence intervals. The (red) star denotes 
the physical value converted to lattice units using = 1588 MeV |69]. We assign an additional 
2% error to this point to estimate the uncertainty in the scale setting method. In Figure |2]^a) 
we display the fit to the lightest two points and in (b) the fit to all four points. Note that the 
extrapolation describes the values of /jr very well. Additionally, one sees that using f^^ or r\ to set 
the scale results in agreement in the extrapolated values, as first observed in Ref. |70j.^ 



^ The scale of r\ is determined through the static-quark potential by solving for r\F(r\) = 1; the values of ri/a can 
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0.1 



0.2 



0.3 



0.4 



0.5 



mj{2y[2nf,) 



FIG. 3: Spin- averaged mass of rjc and J /"^ on the different ensembles. The blue points and purple points 
indicate the masses at mi and m2 respectively. The red line indicates the experimental value. The left panel 
displays the results from the lattice spacing a^^ — 1588 MeV used on all ensembles. This method was used 
to tune the charm-quark mass on the m007 ensemble. The right panel displays the masses scaled by /^r on 
the lattice and extrapolated to f^^^^, as discussed in the text. 

B. Charm-Quark Mass Extrapolation 



To tune the charm-quark mass we use the spin-averaged J/^'-r/c mass. We use the lattice 
spacing determined by MILC {a^^ = 1588 MeV [69j) on the m007 ensemble to estimate the two 
charm-quark masses used for our charm quark propagator calculations. These same two charm 
quark masses, mi and m2, were used on all ensembles. On the MILC ensembles, the value of /3 
was slightly varied for the different light-quark masses. Therefore, the corresponding value of the 
critical mass changes from ensemble to ensemble, leading to a slightly different charm-quark mass 
tuning. This can be clearly seen in the left panel of Fig. [3j where we display the spin-averaged 
J/^'-r/c mass as a function of the light-quark mass, determined with the = 1588 MeV scale 
setting. Ensembles m007 and mOlO share the same value of /3 and therefore the difference in these 
points (the left-most two sets of masses) is due entirely to light-quark contributions, whereas the 
m020 and m030 ensembles each have a different value of (3, so that the variation of the spin-averaged 
mass is due both to light-quark effects as well as a shifted value of the critical mass. 

In the right panel of Fig. [Sj we display our preferred method of determining the charm-quark 
mass using f^^ to set the scale. On each ensemble, we take the spin-averaged J/^-r]c mass and 
divide by the corresponding value of f]^^^ calculated on that ensemble. We then use the value of 



be found in Ref. [71] • 
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li determined in Sec. IV A to scale these values to determine the ratio with f^^"^^, 



It is these scaled values that are plotted in the right panel of Fig.[3]and which we use to extrapolate 
our spectrum calculation to the physical charm-quark mass point, which we take to be 



4/. 



phys 



23.47 , (14) 



phys 



with —!r— = 1.056. (15) 

Jtt 

Here, m^^'^ is taken to be the isospin-averaged pion mass, while f^^'^'^ is taken to be the charged- 
pion decay constant [72]. On each ensemble, we linearly extrapolate the spin-averaged J/^-rjc mass 
(scaled by f^^^^) to the experimental value to determine the parameter mo = rric^^^ (the masses 
of all hadrons are then extrapolated linearly to this charm-quark mass on each ensemble). The 
uncertainties of the extrapolated hadron masses are evaluated using the jackknife method. As a 
check of systematics, we perform the same procedure using the lattice spacing = 1588 MeV 
to perform the linear charm-quark mass extrapolation. Using this second approach, the result- 
ing charmed baryon spectrum is consistent with that of our preferred charm-quark mass-tuning 
method. 

To test the viability of our choice of mixed-action and to gauge the discretization errors, we 
compute both the J/^-rjc hyperfine mass splitting as well as the low-lying charmonium spectrum 
of the XcO) Xci and he- The interpolating fields used for these charmonium states are'^ 

XcO = QcQc, (16) 
xii=QcYl5Qc, (17) 

3 3 

/ij. = ^^e.,fcQc7VQc, (18) 

j=i k=j 



To extrapolate these charmonium masses to the physical light-quark mass values, we use Eq. (11) 
both in quadratic (in mjr) as well quartic form, i.e. 

fn fn ^ 4^ /2 + (4vr)2 /4 ' ^ 



The results of the extrapolation are displayed in Fig. HI and tabulated in Tab. VII In the table. 



^ One can also use improved interpolating operators to extract charmonium states in lattice calculations, especially 
for the excited states Xco, Xci and he, see, for example, Ref. |73| . 
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mJd^flTiQ m,/(2V2";rA) mjd^flnf,) 



FIG. 4: The masses of XcQ, Xci and he as functions of m^/f^. The blue pomts are our numerical values. 
The blue shaded regions show the standard deviation allowed regions of linear fit. The pink shaded regions 
show the standard deviation allowed regions of quadratic fit. The red points are experimental values. 





M^.o (MeV) 


M^^, (MeV) 


Mh^ (MeV) 


Extrapolated Values 


3465(20)(13) 


3525(20)(6) 


3553(25)(14) 


Experimental Values 


3415 


3511 


3526 



TABLE VII: Low-lying charmonium spectrum of XcOi Xci E^nd h^- The experimental values are taken from 
the Particle Data Group [72) . 

the first uncertainty is statistical and the second is an extrapolation systematic from the two 
extrapolation functions used. 

A more stringent test of discretization errors is the calculation of the hyperfine splitting. The 
hyperfine splitting is obtained by fitting the ratio of the two-point correlation functions of J/^ 
and rjc 

to a single exponential 

-jl ^ Ae^^^\ (21) 

where is the mass splitting between the J/^' and r^c- The splittings are first extrapolated to the 
physical charm-quark mass for each ensemble and then extrapolated to the physical light-quark 
mass. As with the charmonium spectrum, we perform a light-quark mass extrapolation using 
both a quadratic and quartic form of Eq. (11). In Fig. [5] we display this extrapolation, finding 



Mj/ij, — M^^ = 93(1)(7) MeV. The first uncertainty is statistical while the second is a systematic 
from the chiral extrapolation. 

It is well known that the lattice computations of the charmonium hyperfine splitting (exper- 
imentally measured to be 117 MeV) are sensitive to the lattice spacing. Qualitatively, one can 
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FIG. 5: Extrapolation of the hyperfine splitting. The blue points are the lattice data. The red point is the 
experimental value. The blue band is the quadratic fit with Eq. while the pink band is the quartic fit 



with Eq. 11 



understand this by performing a Symanzik expansion of the heavy quark action, reveaUng dimen- 
sion five operators arising from discretization effects, which are otherwise identical to the heavy 
quark effective theory (HQET) |74| I75l [76] operator responsible for the hyperfine splitting^ 



(J B 



2 m, 



(22) 



where h^c'^ is the heavy quark field. In the heavy quark action we are using, the coefficients 
of the operators Sb ^ and S'^; Q have been given their tree-level, tadpole improved values 
in order to mitigate the effects of this unwanted discretization effect. It is known the operator 
Sb ^ has a significant effect on the hyperfine splitting [9l [TH I12j . A nonperturbative tuning 
of the coefficient cb can improve the hyperfine splitting in a fixed-lattice spacing calculation; see 
Ref. [77], in particular Fig. 3. However, the qualitative aspects of this effect remain even after 
tuning the coefficients. Previous quenched calculations of the hyperfine splitting have generally 
been low, being about 80 MeV, and showed a strong lattice-spacing dependence. Further, a recent 
direct calculation of the disconnected diagrams has ruled out these (or their lack thereof) being 
the cause of the discrepancy [78 . Our results are consistent with those of the Fermilab/MILC 
Collaboration, which utilized a similar heavy quark action, the same dynamical ensembles and 
staggered light quarks [79]. The Fermilab/MILC Collaboration also performed calculations on 
different lattice spacings, finding similar lattice-spacing dependence to Ref. [79j. Therefore, the 



A proper treatment of heavy quark discretization effects is more involved and can be found in Ref. [9]. 
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discrepancy of our calculated hyperfine splitting with the experimental value is expected. 

C. Light-Quark Mass Extrapolation 

1. Heavy-Hadron xPT Extrapolation 

To perform the light-quark mass extrapolation, we begin with a continuum HHxPT extrapola- 
tion of the baryon masses. The mass formula for these baryons containing a heavy quark was first 
determined in Ref. [63] and later extended to partially quenched theories in Ref. [66]. For dou- 
bly heavy baryons, the was formulated in Ref. [M] and later extended to partially quenched 
theories in Ref. [67^. In this work, we perform SU{2) chiral extrapolations of the baryon masses, 
inspired by Ref. [68^.^ To perform the extrapolations, we treat the J = 1/2 and J = 3/2 baryons 
as degenerate, which is valid at this order in HQET/HHxPT.^ The baryons are grouped into their 
respective SU (2) multiplets allowing for a simultaneous two-flavor chiral extrapolation of all masses 
in related multiplets. This allows us, with only four gauge ensembles, to determine all the relevant 
LECs for a given pair of multiplets in a global fit. The first pair of multiplets contains the Ac and 
Sc baryons. Their SU{2) chiral extrapolation functions are given at next-to-leading order (NLO) 
by 

Mk, ^ Mq 1 cl{^i) ml 6g| .F(m^,AsA,^) 

U fol + Sf{mJU) 47r (47r)2 /i? ' ^ ^ 

Ms, ^ Mo + Ag| 1 c^Mml 

2 g| J^(m^, -AsA,/i) 4 ffl jr(m^,0,/i) 

3(47r)2 /3 +3(4^)2 /3 • ^ ) 

The chiral functions are 



■FK A,^) = (A2-^2^,^)3/.i, [ A + ^/^E^ +E] 3^^.i, f -^)-A3ln ( ^) . (25) 



with 



J^(m,0,^) = W, (26) 



^ For further discussion on SU{2) chiral extrapolations of hadron states with strange valence quarks, see Refs. [801 

MM- 

^ It would be more desirable to use the lattice-calculated masses of the J — 3/2 baryons, but we do not have them 
for this work, and so we use this approximation for now. 
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Fit Range 










5l 


9i 


dof Q 


m007-m030 


1.46(10) 


17.9(2) 


-0.8(5) 


0.2(1.2) 


0.8(1.0) 


-0.1(1) 


0.32 3 0.95 



TABLE VIII: Fit to Ac and Ec masses with NLO continuum formulae. 



< 




(a) 



(b) 



FIG. 6: NLO HHxPT extrapolation of Ma, and Ms, (a) as well as Ms, - Ma, (b). 



and 



:F(m, -A, A) = < 



-J^(m,A,A) +2i7r(A2 



')^/^, m < |A| 



-J^(m,A,A) +27r(m2 - A2)3/2^ m > |A| 



(27) 



To stabilize the fits, we first fit M^^ — to a quadratic in m-,^/ f-,^, and feed this into a fit of the 
masses, yielding the results in Table VIII and extrapolations displayed in Figure |6] One observes 
that the continuum HHxPT fits describe the lattice data very well. However, only the leading 
term, Mq is well determined,'^ while the rest of the LECs, most notably the axial couplings, ^sStt 
and ^sAtt are consistent with zero. This phenomenon is not unique to the charmed baryons. In 
Ref. [53], chiral extrapolations on the nucleon mass in which the nucleon axial coupling, g-i^NN 
(commonly denoted as gA in baryon xPT) was left as a free parameter, returned values which 
were inconsistent with experiment and phenomenology. In fact, given the lattice results for the 
nucleon mass as a function of m^r, it was found that the nucleon scales linearly in m-j^. Such behavior 
signals a delicate cancelation between different orders, a trend which is found in all 2 + 1 dynamical 
lattice computations of the nucleon mass [83]. Therefore, our findings for the axial couplings of 
the charmed baryons are not surprising in this light. To improve the situation, a simultaneous fit 



To determine Mo/f^^" we take our results for Afo//o and scale them by [1 + 5 f [ml^^'' / ff"^"")]-^ . 
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Fit Range 








9l 


dof Q 


m007-m030 


0.85(6) 


19.4(2) 0.6(6) 


1.3(1.2) 5.9(3.9) 


-1.0(6) 


0.04 3 1.00 



TABLE IX: Fit to and ^' masses with NLO continuum formulae. 



[I] 




[I] 16 - 



m^l {2 V2";r f„) m„ / (2 yfln f„) 

(a) (b) 
FIG. 7: NLO HHxPT extrapolation of Ms, and Ms>^ (a) as well as M^^ - Me, (b). 

of the axial charges themselves, along with the masses will most likely be necessary. 

We perform a similar analysis for the J = 1Z2 Hc-H[, isospin doublets, the results of which are 
collected in Table IX and displayed in Figure [Tj The extrapolation formulae for M=/^ and Ms^ 



are similar to those for M^^ and Ma^. They can be deduced by comparing Eqs. (23) and (24) to 
Ref. [66] ,8 

Mh, _ Mo 1 c^fi)ml 3 gi ^(m^, Ah'h, A^) 



U fol + 6f{mJU) 47r 2 (47r)2 /3 



(28) 



Ms' Mo + aSI 



1 4,(^) 

/o l + 5/(mJ/0 47r /2 
1 JP"(m^, -AH'=,/i) 1 5i ^(m^,0,/i) 



2(47r) 



f3 



+ 



2 (4vr)2 /3 



(29) 



The masses of the remaining J = 1Z2 charmed baryons, M=^^, Mq^ and Mq^^, can be treated 
independently. The extrapolation formula for Af=^^ is similar to that of M^,^. There is an axial 
coupling g~.ccB.ccTr well as gE'c^ccT^ where the second coupling is the axial transition coupling of 

* In 5(7(3) HHxPT, the axial couplings for the Hc-H^ system are the same as those for the A^-Ec system, g2 = 
gssn ~ gs's'ir and = gsA-u = gs'Sn- However, in the 5(7(2) theories, they differ by SU{3) breaking corrections. 
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Fit Range 










dof 


Q 


m007-m030 


28.1(2) 


1.4(1.0) 


-1.7(1.0) 


3.0 


1 


0.08 



TABLE X: Fit to J = 1/2 mass with the NLO continuum heavy-hadron formula. 



30 
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FIG. 8: NLO HHxPT extrapolation of M^^^. 

the J = 3/2 to the J = 1/2-tt state. The heavy quark symmetry also requires these couplings to be 
the same in the heavy quark limit. At this order, we can treat the J = 3/2 E*^ as degenerate with 
the Sec- The results are collected in Table [X] and displayed in Figure [Sj with the extrapolation 
formula j67] 

where we have set A=*= = in this analysis, valid at this order in the heavy-quark expansion. 

One feature which is more pronounced in this fit is g^ < 0. Taken at face value, this would 
suggest the Lagrangian was non-Hermitian, and the theory not sensible. Therefore, even though 
these fits reproduce the lattice data well and predict a mass within a few percent of the physical 
value, they must be taken with caution. Most likely, as with the nucleon mass [83j, there is a 
delicate cancelation of terms at different orders, and therefore one does not have confidence in 
these determinations of the LECs. 

Similar to the s = —3 $7, the J = 1/2 0,c and Qcc do not have mass corrections which scale as 
m^. This is because these baryons do not contain any valence up or down quarks, and therefore, the 
leading SU{2) axial coupling vanishes |681 18^ . The SU{2) chiral extrapolation formula for these 
baryon masses is then expected to be as convergent as that for pions. The mass extrapolation 
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n 


Fit Range 






(4) 

"a 




x' 


dof 


Q 




m007-m030 
m007-m030 


20.4(6) 
27.7(4) 


-3.0(4.6) 
-7.3(3.0) 


46(61) 
109(40) 


-164(227) 
-392(149) 


0.00 
0.00 









TABLE XI: Fit to J = 1/2 ilc and ilcc masses with NLO continuum heavy-hadron formulae. 




u 26 - 



0.5 0.0 




0.2 0.3 



(a) (b) 
FIG. 9: NLO HHxPT extrapolation of Mn, (a) and Mn^^ (b). 



0.4 0.5 



formula for the i7c and ilcc are both given by 



Mn Mo 



1 



c- m: 



- + 



mz, 



+ ^(4)(.) 



(31) 



At this order, the two-loop corrections from should be included as corrections to and /3q • 
Further, there is a ln^(m7r) correction with fixed coefficient. However, since we only have four mass 
points, we cannot judge the quality of the fit anyway, so we ignore these corrections. The results 
are collected in Table XI and displayed in Figure |9j Performing a fit with an = and (3n = 
returns consistent mass predictions with smaller uncertainties. We take the zero-degree-of-freedom 
fit as our central result as it provides a more conservative uncertainty. 



2. Polynomial Extrapolation 



Given the issues of performing the heavy-hadron chiral extrapolations as discussed above, we 
also perform polynomial extrapolations in m^. We use the difference between the polynomial 
extrapolations and the heavy-hadron chiral extrapolations as an additional estimate of systematic 
extrapolation uncertainty. We use up to three different polynomial fit functions for each of the 



21 





1 08 

1 .v/O 




i .WO 


Cl 

X 


1.04 


s 






1.02 


^~ 


1.00 




0.98 




0.96 



Ac 
s = 
c=l 



s = 
c= 1 



i =-1 

c= 1 



s=-l 
c=l 



s = 
c = 2 



He 
s=-2 
c= 1 



FIG. 10: Ratio of extrapolated masses to experimentally measured masses. The first point represents the 



HHxPT fit, the second point is a fit with Eq. (32 1, the third with Eq. (33 1 and the fourth with Eq. (34) 



State 
(J =1/2) 



Ma, 
[MeV] 



Me, 
[MeV] 



Me. 
[MeV] 



Ms^ 
[MeV] 



Ma, 
[MeV] 



Mh„ 
[MeV] 



Ma„ 
[MeV] 



Prediction 2342(22)(11) 2527(17)(13) 2527(20)(08) 2638(17)(10) 2687(46)(16) 3665(17)(14) 3680(31)(38) 
Exp. Mass 2286 2468 2454 2576 2698 3519 



TABLE XII: Direct light/heavy quark mass extrapolation of the J = 1/2 charmed baryon spectrum. 



charmed hadron masses: 



M2 _ Mq 1 

U /o 1 + Sfim^/fTr) 

M3 Mo 1 
X 

M4 _ Mq 1 



/o 1 + SfimJU) 



2 
J n 










ml 

J TT 




73" 


+ C3 


ml 

J TT 




4 


+ C4 



(32) 
(33) 
(34) 



/o 1 + 5f{mJU) 

In Figure 10 we display the results of these fits as well the heavy-hadron fits as ratios with 
respect to the experimental masses. The experimental values for the baryon masses are taken from 
the Particle Data Group [72 . As it can be seen, there is very little variation in the results of the 
extrapolated masses. In all cases, the different extrapolations are consistent within one sigma. 



In Table XII we provide the extrapolated baryon masses, taking the central value from the 
HHxPT extrapolations. The first uncertainty is statistical and the second uncertainty is a com- 
prehensive systematic uncertainty. This systematic uncertainty is derived by comparing the poly- 
nomial light quark mass extrapolations to the HHxPT extrapolation. Further, it includes the 
uncertainty associated with the choice of fitting window for the correlators as well. Except for the 
Oc, the extrapolated masses are systematically high, indicative of a discretization error. 
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D. Discretization Errors and Mass Splittings 



In this work, we have performed calculations at only a single value of the lattice spacing, 
with a ^ 0.125 fm, prohibiting us from performing a continuum extrapolation. However, we can 
take advantage of various symmetries and power counting to make a reasonable estimate of the 
discretization errors present in our calculation.^ In these heavy-light systems, the discretization 
errors arise both from the light and heavy quark actions. The corrections from both generically 
scale as 0(0^) for each of the charmed baryon masses. If we consider SU{3) symmetry, then the 
leading discretization errors for all baryons in a given SU(3) multiplet must be the same, with 
corrections scaling as 0(0^ {nis — mu))- Further, if one considers the combined large- A'^c, SU{?>) 
and heavy-quark symmetries |85J, then all the singly charmed baryon masses we calculate in this 
work share a common discretization correction to their masses, with sub-leading corrections scaling 
as 0{a^ /Nc) as well as the SU{?>) breaking corrections. Therefore, all the singly charmed baryon 
masses we compute in this work, {Ac, He, So ^c} share a common discretization correction, 
which happens to be the dominant discretization error. The same analysis holds for the doubly 
charmed baryons as well, {Hcc, i^cc} with a common error, albeit different from the singly charmed 
correction. It is therefore advantageous to consider extrapolations of baryon mass splittings, as 
these mass splittings exactly cancel the leading discretization errors. 

Before proceeding with the analysis of the mass splittings, we first use power counting arguments 
to estimate the discretization errors. The leading discretization corrections from the light and heavy 
quark actions can be estimated as [13] 



where p is a typical momentum scale, of the order of Aqcd, the characteristic hadronic scale. To 
be conservative, we can take Aqcd = 700 MeV which leads to the estimates 



^ With a single lattice spacing, we can not disentangle both the discretization errors and the tuning of the charm 
quark mass. The effects we discuss here as discretization errors are really a combination of the two. 

" With the full J — 3/2 and J — 1/2 heavy baryon mass spectrum, one could perform an analysis of the large- A^c 
baryon mass relations [86J W7\ as has recently been performed for the light quark octet and decuplet baryons [88] . 




(35) 



5g(a2) = 68 MeV, 
dqia^) = 19 MeV. 



(36) 
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FIG. 11: Polynomial extrapolations of of J = 1/2 mass splittings amongst lieavy-quark-S'f7(3) multiplets 



with Eq. (341 



When considering mass splittings amongst a given 5f7(3) multiplet, these leading errors become 
further suppressed by — mu effects, 



AA M (n^\ - ^s{nic)iap) ■ ml -ml 
) - 2(1 + am,) ^'^'^^^ A| 



(37) 



Mass splittings between the two singly charmed SU{3) multiplets, AM^'^, would receive simi- 
lar discretization corrections, with the extra suppression of l/Nc- Combining these estimates in 
quadrature, we estimate the discretization errors for the baryon masses, and various mass split- 
tings (using A;^ = 2\/27r/7r and the physical kaon and pion masses) 

5Mh^ = 71 MeV , 
5Mh^^ = 78 MeV , 
5AMh^ = 12 MeV , 
5AMh,^ = 13 MeV , 

7-6,3 



dAM^l" = 24 MeV , 
SAM^'^ = 26 MeV . 

'f'CC 



(38) 



Given our limited number of light-quark mass values, we are not able to perform the (mixed- 
action) HH^PT analysis of the mass splittings. We therefore perform our fits using the polynomial 



For the doubly charmed baryon masses, we double the estimated heavy quark discretization error. As mentioned 
above, this uncertainty also includes any miss-tuning of the charm quark mass, and thus a double charmed baryon 
will be miss-tuned twice as much. 
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FIG. 12: Ratio of extrapolated mass splittings to experiment [72]. The first point is a fit with Eq. (32), the 



second with Eq. ( 33 1 and the third with Eq. ( 34 1 



fit functions, Eqs. (32)-(34), with Mq replaced by Aj|^^^^. We perform the extrapolations of the 
mass splittings, Ms, - Ma,, {Ms^,MnJ - Ms,, Ms, - Ma, and Mq„ - Mh„. In Figure 



11 



we 



display the extrapolation of these mass splittings using Eq. (34) and in Figure 12 we show the ratio 



of these fits to the experimental values. Our final predicted splittings are determined by using the 
quartic fit function as the central value with the differences from the quadratic and cubic fits to 
estimate light quark mass extrapolation errors (in addition to those from the quartic fit). 

As discussed earlier in this section, the dominant discretization error in the mass calculations 
is common to all baryons, given the various symmetries. Therefore, this correction will shift all 
the baryon masses in one direction. We can determine the sign of this correction in the following 
manner. First, we can determine the singly charmed baryon spectrum by taking our extrapolated 



mass splittings, column (a) of Table XIII, and using M^J^ and M^J^ as reference scales, M^^^^ = 
M^^^ + AM/i^__A^.s.i resulting in the predicted masses. Table XIII (h). We then compare these to 



^A, s, ^^''-'/ic~A,,s,) resulting in the predicted masses. Table 
our direct mass extrapolations M^^'''^'^*, given in Table 



XII 



The first method is free of the leading 



discretization errors while the second is not. We can then construct the quantity. 



fAivcct 



(39) 



which is a measure of these discretization errors. The sum runs over all four singly charmed baryons 
he for which we have both methods to determine the masses (A'^^, = 4). The first thing to note is 
that every element contributing to the sum is a positive quantity, suggesting the discretization errors 
increase the baryon masses. It is also interesting to note that in our calculation, 6Mc{a^) = 59 MeV, 



comparable to our estimated leading discretization effects, Eq. ( 38 ) . We can then refine our estimate 
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State 


Latt. Pred. 


Exp. 


State 


Mass Split. 


Direct Mass 


Exp. Mass 




[MeV] 


[MeV] 




[MeV] 


[MeV] 


[MeV] 








Ma. 




2342±22± lli?i 


2286 


Mh. - Ma. 


164 ±14 ±23 ±12 


182 


Mh. 


2450 ± 14 ± 23 ± 12 


2527±17±13t?i 


2468 


Ms. - Ma. 


190 ± 27 ± 18 ±27 


168 


Ms. 


2476 ± 27 ± 18 ± 27 


2527±20±8t?i 


2454 


Mh- - Ms. 


113± 18±8± 12 


122 


Mh^ 


2567 ± 18 ± 8 ± 12 


2638± 17± 10i?i 


2576 


Mq^ - Ms. 


195±21±7± 12 


244 


Mn. 


2649 ± 21 ± 7 ± 12 


2687±46± 16i?i 


2698 








Mh.. 




3665± 17± 14t?g 


3519 




98 ± 9 ± 22 ± 13 




Mn^^ 


3763 ± 19 ± 26 tjl 


3680±31±38i?8 





(a) 



(c) 



TABLE XIII: Resulting charmed spectrum, extrapolated in the light-quark mass to the physical m^^y^ / f^^^^ 
point. In (a) we display the mass splittings of the baryons related by SU{3) and large Nc symmetry. As 
discussed in detail in the text, the first uncertainty is statistical, the second is systematic and the third 
is our estimate of discretization errors. These are the central results of this work. In (6), we display our 
resulting baryon spectrum determined using the experimental values of M"^^^ and , combined with our 
splittings in (a). For the flee, we use our extrapolated value of Afs^^ given the present uncertainty in the 
experimental value. In (c), we present the results of our direct mass extrapolations, including our estimated 
discretization errors. The results from the two methods are consistent at the one-sigma level. 

of the leading discretization errors to be 



tn MeV , 



(40) 



where we have also assumed that the doubly charmed discretization errors do not change sign 



relative to the singly charmed baryon corrections. Our final numbers, collected in Table XIII 
include these discretization error estimates in the quoted uncertainties. 



DISCUSSION AND CONCLUSIONS 



The central results of this work are the predicted mass splittings, displayed in the left panel 



of Table XIIT The first uncertainty is statistical and the second uncertainty is a comprehen- 
sive systematic as discussed in the text. The third uncertainty is an estimate of discretiza- 
tion errors, which must scale as 0{a'^{ms — rriu)) for members of the same SU{?>) multiplet or 
0[a? /Nc) + 0{a^{ms — m„)) otherwise, as dictated by the approximate symmetries. These results 
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FIG. 13: Comparison among charmed baryon mass splittings of dynamical lattice calculations. The results 
of Na et al. are taken from Ref. [25]. 

have been extrapolated to the physical charm quark mass and the physical light quark mass defined 
respectively by 

r — = 23.47 , 

4jphys 

phys 

^ = 1.056. 

Jn 

To perform these extrapolations, we first formed the dimensionless ratios {Mj^^^ — Mj^^^) / f]^^^ , 
taking into account the known light-quark mass dependence of /jr. The mass splittings in MeV 
are then determined with /^r = 130.7 MeV. These physical values are all taken from the PDG |72] . 



In Fig. 13, we compare some of our mass splitting results with those of Gottlieb and Na 
the only other dynamical calculation of the charmed baryon spectrum. They used the same MILC 
gauge ensembles, as well as the fine a ~ 0.09 fm lattices. For the light quark propagators, they 
used staggered fermions, and for the heavy quark, an interpretation of the Fermilab action was 
used, defining the charm mass with the kinetic mass instead of the rest mass. Their work is still 
somewhat preliminary and does not yet provide a systematic uncertainty. However, our results are 
consistent with theirs, especially those on the same ensembles with a ~ 0.125 fm. 

We additionally use these mass splittings, combined with the experimental value of M"^^ and 
M^^P to determine the J = 1/2 baryon masses. Aside from the Hcc state, the masses determined 



Because the Ecc has not been verified by muhiple experimental groups O |6l [721 [89l [90], we chose to use our 
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FIG. 14: A summary of charmed baryon masses in MeV calculated using LQCD. We show both of our 
methods for obtaining the spectrum, the direct mass extrapolation (Liu et al. 1) and also using the extrapo- 
lated mass splittings, combined with M'^^ and AI^^ (Liu et al. 2). These results are taken from Table 
The other results, displayed for comparison, are taken from Table pOV] 



XIII 



in this way are consistent with our direct mass extrapolation results, Table XIII (c), after including 
our estimated discretization errors. We used power counting arguments |13| l9T] to estimate the 
size of these corrections and we compared our two methods of determining the baryon masses to 
determine the expected sign of the leading discretization corrections. In Fig. [T4j we display our 
resulting mass calculations using the results from both the mass splitting method (Liu et al. 2) as 
well as the direct extrapolation of the masses (Liu et al. 1). Additionally, we compare these with 



results from previous calculations, found in the Refs. of Table XIV (for those calculations with 
more than one lattice spacing, we show only the results from the ensemble with lattice spacing 
closest to the one used in this work). 

Finally, we compare the doubly charmed baryons with the predictions of theoretical models, as 



shown in Fig. 15 Although the SELEX Collaboration has reported the first observation of doubly 
charmed baryons, searches by the BaBar [89], Belle |90] and Focus [95j Collaborations have not 
confirmed their results. This makes it interesting to look back to the theory to see where the 
various predictions lie. We compare with a selection of other theoretical results, such as a recent 
quark-model calculation [96], relativistic three-quark model [97], the relativistic quark model [98] . 



extrapolated value of Ms^^, combined with our extrapolated value of A/n^c ~ ^^2cc to make a prediction for the 
Qcc mass. 
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Group 


Nf 


S'h 


ar^ (GeV) 


L (fm) 


Bowler et al. [19] 





tree clover [92| 


2.9 


1.63 


Lewis et al. [22j 





D234 [§3] 


1.8, 2.2, 2.6 


1.97 


Mathur et al. [23] 





NRQCD [92] 


1.8, 2.2 


2.64,2.1 


Flynn et al. [20] 





NP clover 


2.6 


1.82 


Chiu et al.[2T] 





ODWF [94] 


2.23 


1.77 


Na et al. piling] 


2 + 1 


Fermilab [5] 


2.2, 1.6, 1.3 


2.5 


This work 


2 + 1 


Fermilab 


1.6 


2.5 



TABLE XIV: Summary of existing charmed baryon published calculations from lattice QCD. Please refer 
to the above references and references within for more details. 

the heavy quark effective theory [99] and the Feynman-Hellmann theorem [100| . We compute the 
mass of Hcc to be 3665 it 17 it 14^73 MeV, which is higher than what SELEX observed, although 
less than two sigma with our estimated discretization errors; most theoretical results suggest that 
the Hcc that is about 100-200 MeV higher than the SELEX experimental value. To improve this 
situation, we need results at multiple lattice spacings to reduce this systematic uncertainty. The 
Occ mass prediction made by this work is 3763 it 19 it 26 MeV, and the overall theoretical 
expectation is for the 0,cc to be 3650-3850 MeV. We hope that upcoming experiments will be able 
to resolve these mysteries of doubly charmed baryons. 

Our largest uncertainty presently arises from the lack of a continuum extrapolation. Therefore, 
in the future we plan to extend these calculations to a second lattice spacing. This will hopefully 
allow us to significantly reduce the size of our discretization errors. Additionally, we are extending 
our calculation to include the spin-3/2 spectroscopy. 
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